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I. INTRODUCTION 



The problem of quantitative description of multicomponent systems is interesting from 
both theoretical and practical points of view. The theoretical description of real systems from 
"first principles" is now impossible because of absence of constructive methods of partition 
functions calculations for realistic models. Therefore the only way of the realistic models 
investigations is the phenomenological approach based on some additional assumptions. 
One of the most popular phenomenological approaches to statistical thermodynamics of 
condensed systems is the Cahn-Hilliard and Ginzburg-Landau approximation. Application 
of these very attractive models has some well known restrictions. Most essential of these 
restrictions is absence of concrete relations between models parameters and characteristics 
of constituents and their interactions. 

At construction of the theory of these structures there is a lot of the difficulties caused by 
absence of constructive methods of the analysis of realistic models of substance from "first 
principles". Therefore it is necessary to be content with development of phenomenological 
models. The examples of such models are Ginzburg-Landau tl|, Cahn-Hilliard y, 0] 
models and the various variants of lattice and cell models [sl, H, 0, @] . 

The basic idea of lattice models consists on introduction of some ideal lattice on which 
sites atoms of components are distributed. This idea allows to consider the identical short- 
range repulsions between atoms, but does not allow to consider differences in repulsions of 
the components. Most essential disadvantages of the lattice models are well known: 

1. The assumption of some lattice existing independently of difference between compo- 
nents proper atomic sizes. This assumption is comparatively plausible by the condition 
of vanishingly small difference of the components atomic sizes. As a rule, this condi- 
tion takes no place. Really, introduction of lattice is the way the short-range parts 
of interatomic potentials account. Difference of short-range parts of components in- 
teratomic potentials leads to lattice distortions and the notion of the lattice, strictly 
speaking, in this case becomes invalid. 

2. The postulate on Fermi-like functional form of distribution function for average oc- 
cupations numbers has not any theoretical justification — there are infinitely many 
functions with range of values [0; 1]. 

3. Interatomic potentials in lattice-like models of condensed systems, as a rule, cannot 
be described with small number of parameters, such as interaction energy of nearest 
neighbors interactions. The unavoidable lattice distortions lead to the changes of 
interatomic distances and short-range parts of interatomic potentials are not slowly 
changing functions. 

4. The assumption on pair-wise interatomic potentials is not very convincing, because 
all the interatomic interactions appear as result of the averaging over fast electronic 
degrees of freedom in system. These degrees of freedom are very sensitive to atom 
coordination environment. 

Ginzburg-Landau and Cahn-Hilliard models contain a set of phenomenological parame- 
ters which connections with real physical characteristics of the components and their inter- 
actions remains hidden. 

The generalized lattice model (GLM) of multicom pon ent condensed systems (such as 
solid or liquid solutions) was proposed in papers B HEl and developed in |12l . Il3l . |14| . In 
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contradistinction to usual lattice model (see for example jl, 0, 0, S]), the GLM takes into 
account the following factors: 



The short-range interatomic repulsions |9|]. These repulsions are not identical for 
different pair of atoms, therefore it is impossible to take into account the repulsions by 
means of lattice introduction (at best, the lattice model can be considered as method 
of short-range interactions account in the case of one-component system and that with 
some essential restriction, related in particular to thermal defects). 

The presence of the local fields due to the long-range parts of the interatomic poten- 
tials These fields have the essential influence on both equilibrium properties [§] 



and non-equilibrium processes [ll(| on the corresponding scales. 



3. The connection between the GLM and the Ginzburg-Landau approximation is ob- 
tained in paper The relation of some characteristics of the components and 
their interatomic potentials with Ginzburg-Landau parameters is established. This 
relation makes possible to use the mathematical tools of the Ginzburg-Landau and 
Cahn-Hilliard approximations for the GLM research. 

4. The existence of comparatively stable polyatomic complexes that manifest itself in 
both thermodynamics and kinetics as one and indivisible particles [l3| . 

The present paper contains further development of the GLM. In addition to the previous 
results it takes into account the following factors: 

1. Existence of the internal atomic degrees of freedom such as atomic electric and mag- 
netic moments. These degrees of freedom are responsible for the local magnetization 
and local electric polarization in the system. 

2. Reduction of the GLM to the Ginzburg-Landau-Cahn-Hilliard-like (GLCH) functional 
form. In contrast to the GLCH theory, this functional contains the well defined pa- 
rameters that have direct connections with the characteristics of the constituents and 
their interactions. 

3. Colossal times of the structure transformation and the equilibrium reaching. Real 
condensed systems rule, essentially non-equilibrium systems. 

4. The presence of the hierarchy of relaxation times in real condensed systems and related 
partially equilibrium states of the systems. The evolutional processes rearrangements 
in such systems have multi-stepped character. 

5. Absence of full thermodynamic equilibrium in the system. 

All the results are formulated on the basis of the unified mathematical apparatus and the 
common physical ideas. All the approximations have clear physical sense, well foundations, 
and strictly based conditions of the apphcability. 
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II. GENERALIZED LATTICE MODEL THE BASIC NOTIONS AND RELA- 
TIONS 



A. Conditions for free energy 

The basic idea of the lattice models is the assumption that particles are located in sites 
of some ideal lattice. This assumption is incompatible with differences of atomic sizes of the 
components and with presence of the various kinds defects in real condensed systems. 

Therefore let us introduce instead of the ideal lattice the connections between local den- 
sities of the components with a view to take into account the differences of the atomic sizes 
of the components and the various kinds of defects in the system. 

It should be noted that the short-range repulsive parts of the interatomic potentials lead, 
in particular, to some restriction on the local densities nj(r) of the constituents particles 
numbers in the system 

n^{r) < -, (1) 

where uJi is the inverse value of the maximal local density of z-th component (i = 1 -i- m, 
m is the number of the components in the system). The quantity uji has dimensionality of 
volume and henceforth will called as the specific atomic volume of i-th component. As far 
as the quantity Uini{r) is the local volume fraction of i-th component at the point r, then 
we have the restrictions on the local densities of the components for all points in the system: 



J2 Mr) < 1- (2) 



1=1 



Let us introduce the additional component — the vacancies (or holes) with their proper 
volume Uq and local density no(r). Suppose the holes do not interact with real components 
but they fill all the unavoidably existent vacant places in the system (such as thermal or 
radiation defects). With account of the vacancies we have the following the packing condition: 

m 

ni{v) -1 = 0. (3) 

1=0 

The minimization of free energy should be realized under the packing condition, which takes 
into account the short-range parts of the interatomic potentials in the system. Hence the 
interatomic potentials should be included into the free energy with cutting out their the 
short-range parts: 

K (r\- I 1^1 - ^'^^ (A) 

'^^^ \ 0, otherwise, 

where Wij{r) is "true" interaction potential between i-th and j-th components, Oij are the 
cutting parameters, related to the specific atomic volumes of the components by the relations 



[{u.f' + (CO,)'/'] . (5) 

Beyond the packing condition the numbers Ni of components atoms in the system 
should be fixed at the free energy minimization. These conditions have the following form: 

JnM,.-N.^O, »^0.». (6) 
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The packing condition ([3]) and condition of numbers particles conservation ([6]) should be 
satisfied for any form of the configuration part of the free energy. 



B. Configuration part of free energy 

Denote by Kij (r — r') the long-range part of the independent of the inner degrees of 
freedom (such as their electric and magnetic moments) potential energy of i-th and j-th 
interacting particles located in points r and r', respectively. 

The number of z-th kind particles in an infinitesimal volume dr near a point r is nj(r) dr. 
Similarly, the number of j-th kind particles in infinitesimal volume dr' near a point r' is 
nj{r') dr'. Then the full energy of interactions between these particles near points r and r' 
is 



Hence the inner degrees of freedom independent part energy of the system have a form: 



The similar way leads to expression for configuration part of the free energy with account 
of the atomic electric and magnetic moments of the components: 



where Dj(r) and Mj(r) are electric and magnetic atomic moments of i-th components, 
located at the point r, Kij{r — r'), Qij{r — r'), Rij{r — r') and Sij{r — r') are long-range parts 
of the relevant two-body interatomic potentials. 

It should be noted that the vacancies do not contribute into the configuration free energy, 
but their contribution into the free energy due to entropy term: 



ni{r) dr Kij{r — r') nj{r') dr' . 



(7) 




(8) 




(9) 




m 



njjr) 
n{r) 



) 



S = - ^ ni(r)ln 



dr 



(10) 



where T2,(r) is the summarized local density of the particles and vacancies 



m 



n 



(r) = ^nj(r). 



(11) 



3=0 
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Thus, the Helmholtz free energy F of the system with account of the non-uniform external 
electric E(r) and magnetic H(r) fields has the following form: 



^ = // - r')^.(r)n,(r')drdr'+ 

II ^^^(^ - r')^.(r)^,(r') (D,(r) ■ D,(r')) drdr'+ 
+ i E // ^^^(^ - r')^.(r)%(r') (M,(r) ■ M,(r')) drdr'+ 
+ ^ E // ^^^'•(^ - r')n,(r)n,(r') (D,(r) ■ M,(r')) dTdv'+ 
+ E / (E(r)-D,(r))n,,(r)rfr+5^ / (H(r) • M,(r)) n,(r) rfr+ 

+^E / 



Uiir 



(12) 



n r 



where T is absolute temperature in energetic units. 

The equilibrium space distributions of components is determined by minimum of 
Helmholtz free energy functional with account of the conditions (IHl) and ([6]). But it should 
be noted that beyond this condition we need some additional conditions connecting all the 
electric Dj(r) and magnetic Mj(r) moments with the external fields. Suppose that exter- 
nal fields influence on the orientations of related moments but do not influence on their 
magnitudes. Hence we have 



{Bdr)f-D^ = 0- (M,(r))2-M2 = 0. 



(13) 



C. The Lagrange functional and full equilibrium in the system 

For the minimization of the Helmholtz free energy (fT2!) at the conditions ([3]), ([6]), (fT3|) 
let us introduce the Lagrange functional £ depending on the local densities ni(r) of the 
components, their electric and magnetic moments Dj(r), Mj(r), the external fields E(r), 
H(r), and the Lagrange multipliers ^(r), Aj(r), //^(r), fif 

£({n.(r)}, {D, (r)} , {M, (r)} , {A. (r)}, {u, (r)}, {^{r)}, ^,) = 



m 
i=Q 



nJr) dr - Ni 



(V) 



i=l 



Xi (r) 



[(D, (r))^ - A'] dr - 



i=l 



J^iir] 



(V) 



[{Mi {r)f - M^] dr- f ^(r) (f^ " 1 ) ^^r. 



(14) 
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The necessary condition or the extremum of the Lagrange functional is vanishing of the 
ordinary and functional derivatives with respect to ni(r), Di(r), Mj(r), ^'(r), \i{r), ^'^(r), 
IJ'i- 

< O; 



5ni (r) 
5C 

SMi (r) 
6C 

dC 

dfii (r) 
6C 



Sui (r) 



= 0; 
= 0; 
:0; 
0; 
0; 



0. 



Simple calculations lead to following system of equations: 



Sm (r] 



Pi, = -a;,*(r)+rin(^^) J K,,(r - r')n,(r')rfr'+ 

(V) 

m p 

m „ 

+ i?,,(r-r')(M,(r)-M,(r'))n,(r')(ir'+ 

7=1 //.s 



r — r 



, r(D,(r).M,(rO) + (M,(r)-D,(r'))} 



(15) 



(16) 



nj(r')dr'+ 



+ (E(r)-D, (r)) + (H(r)-M, (r)) ; 



5C 



^D, (r) 

m „ 

y (r) g,, (r - r') n, (r') D, (r') dv'+ 



{V) 



(17) 



f ^^(^) -^^^^ - "^ (^O - ^i(r)Di(r) + E(r)n,(r) = 0; 
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SM, (r) 

m „ 

J2 j <v) R,^ (r - r') n, (r') (r') dr'+ 



(18) 

1 "' f 







5^1 i (r) 
uJi Uiiv) -1 = 0; 

i=0 



(V) 

a, fr) (21) 
(D.(r))^-A' = 0; 





5vi (r) (22) 
[M, (r))^ - = 0. 



This system of integral equations fll6tl22p describes the space distributions of the compo- 
nents as well electric and magnetic moments in the system with account of both short-range 
and long-range parts on interatomic potentials in presence of the external fields. This system 
of equations take place for the case of full thermodynamic equilibrium in the systems. 

Unfortunately, at present there are no effective methods of such kind equations solutions 
for kernel Kij{r), Qij{r), Rij{r), S'jj(r) of general form. But under the some conditions this 
system of equations can be reduced to a system of the partial differential equations. Instead 
of the interatomic potentials Kij{r), Qij{r), Rij{r), Sij{r), the system of partial differential 
equations contains a set of integral characteristics of the potentials. 

III. GINZBURG-LANDAU-CAHN-HILLIARD-LIKE APPROXIMATION 

A. Reduction of the Generalized lattice model to Ginzburg-Landau-Cahn- 
Hilliard-like approximation 

There are at least three scales of the sizes in the system: 

1. atomic sizes Oq; 

2. range of actions of long-range parts of the interatomic potentials ro; 
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3. distances 60 on which changes local compositions and/or local moments of the com- 
ponents in the system. 

Suppose these parameters obey the inequalities: 

ao<ro<^bo. (23) 

Let us transform the expression for moments-independent part of the configuration energy 
in ([12]): 

Pi = lf2 J J - drdr'. (24) 

With regard to condition 

|r-r'|<feo (25) 

we have 

3 n / N 1 3 



n.Ar') ^ n,ir) + ^{x'^ - x.) + ^ - " ^.J, (26) 

s=l 



Sl,S2 = l ^ ^ 



where Xg and x'^ denote Cartesian components of the vectors r and r', respectively. Substi- 
tution of fl26l) into (1211) lead to the result: 



<"^= / K,,ir')\rrdr'. (28) 



where A is the Laplace operator, K^j^ are some integral characteristics of the interatomic 
potentials 

iV) 

The expressions for moment-depending parts of configuration energy in (|T2|) can be trans- 
formed by analogy: 

F,{{n,{r)} , {D,(r)}) = J l^f^ (D.(r)n.(r) ■ A [D,(r)n,(r)]) + 

{V) ^ ''^=' (29) 

+ 2 E (D.(r) ■D,(r))n,(r)n,(r) dr, 
ij=i J 



F3(K(r)},{M,(r)}) = J j^^ E (M,(r)n,(r) ■ A [M,(r)n,(r)]) + 

(V) *'^=^ (30) 

m 1 

+- J2 ^ (M.(r) ■ M,(r)) n,(r) n,{v) dr, 

ij=l ) 



F,{{n,{r)} , {D,(r)} , {M,(r)}) = / ^ (D.(r)n,(r) • A [M,(r)n,(r)]) + 

(V) ^ *'^=^ 



+ i J2 (D,(r) . M,(r)) n,(r) n,(r) I dr, 



(31) 



where -Rjj'* and S^f are the integral characteristics of the moment- depending parts of 
the interatomic potentials determined by perfect analogy with fl28l) . 



Substitute (|28|), ([29]), (|30|) and (|3T]) into ([T2|): 
F(K(r)},{D,(r)},{M,(r)}) = 



iV) 



^ m 

-y 

12 ^ 



+gSf (D,(r)n,(r) ■ A [D,(r)n,(r)]) + (M,(r)n,(r) ■ A [M,(r)n,(r)]) + 



(D,(r)n,(r)-A[M,.(r)n,(r)]) 



+ 



4e 



a1 ' + Ql ' (D.(r) . D,(r)) + ii< ' (Mi(r) • M,(r)) + 



i,j = l 1- 



(D.(r)-M,(r)) 



(32) 



[E(r) ■ Di(r)) n,(r) + (H(r) ■ Mi(r)) ^^(r) + Tni(r) In 



nAr] 



n{r) 



dr. 



The next step is using the Green formula 



u (r) Av (r) dr = — J (Vm (r) ■ Vv (r)) dr, 

(V) (V) 



(33) 



that has place if the functions u{r), v (r) and their gradients Vm (r), Vf (r) vanish on the 
boundary of domain V. Using this formula, the terms with Laplacians in (!32l) take the 
following forms: 



rii (r) Arij (r) dr = — J (Vrii (r) ■ Vuj (r)) dr, 



(34) 



I (D, (r) (r) ■ A [D, (r) n, (r)]) dr = - ^ ^ (V [Dr(r)n, (r)] ■ V [DJ(r)n, (r)] ) rf: 

(V) (V) 



(35) 
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{V) 



(M, (r) (r) ■ A [M, (r) n, (r)]) rfr = - / (V [Mf (r)n, (r)] ■ V [M^{r)n, (r)]) dr, 

(V) 



(36) 



3 

(Di (r) (r) ■ A [M, (r) Uj (r)]) = - / (V (r)] ■ V [M;(r)n, (r)] ) dr, 

(F) iV) 

(37) 

where i^f (r) and M"(r) denote components of the vectors Dj (r) and Mj (r), respectively. 
As a result, we have the following expression for the Helmholtz free energy functional: 



F(K(r)},{D,(r)},{M,(r)}) 



(V) 
3 



-y 

12 ^ 



(Vn,(r)-Vn,(r)) + 



+Qf E (V [A"(r)n.(r)] ■ V [l^-(r)n,(r)]) + 



a=l 
3 



+4'^ E [^r(r)n.(r)] . V [M;(r)n,(r)]) + 



a=l 
3 



E (V [^r(r)n.(r)] • V [M;(r)n,(r)]) 



a=l 



+ 



m 

4e 



a1»' + Q,'°' (Di(r) . D,(r)) + (M.(r) • M,(r)) + 



i,j=l >- 



(38) 



+4°) (D.(r)-M,(r)) 



i=l 



[E(r) ■ Di(r)) ni(r) + (H(r) ■ M,(r)) ni(r) + rn,(r) In 



n,' r 



r?,(r) 



dr. 



This functional is similar to Ginzburg-Landau and Cahn-Hilliard functionals, but in con- 
trast to these functionals formula (l38l) 

1. is not restricted by the polynomial over order parameters in the integrand; 

2. all the parameters in fl38l) have clear physical sense. 

The functional fl38l) has some essential advantages over the more exact functional f|T2l) : 

• The expression f l55]) does not contain any unknown functions of general form (such 
as interatomic potentials), but contains the finite set of their simple characteristics 



(numerical parameters) K,^ 



(p) 
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Analysis and solution of differential equations, that can be obtained from functional 
(!38|) . are much more simple problems as solution of non-linear integral equations (fT6] - 
22|) for general case; 



• The inverse problem — search of parameters using some experimental data (such as 
phase diagrams) — is the immense problem from the integral equations f|T6] - [22l) . but 
quite realistic one from the functional (1381) : 

• At last, the solution, based on the functional fl551) . can be used as well initial approx- 
imation for much more complicated problems, related to functional ( (fT6] - [22|) ). 

There is the necessary condition of the mathematical correctness for the functional 



The necessity of the fluctuations suppression in the system means that all the matrices Kl- , 
Qfji ^if y ^if Klust be negative-definite. 



B. The Lagrange functional for GLCH-like approximation and full equilibrium in 
the system 

It should be noted the conditions ([3]), ([6]), (fT3|) must satisfied not only for functional (fT2l) . 
but also for functional (1381) . Let us introduce the Lagrange functional £i, corresponding to 
the GLCH-like Helmholtz free energy functional fl38p . by analogy with (1141) : 

£i({n,(r)}, {M,(r)} , {D,(r)} , {A^ (r)}, {z/^ (r)}, {^{r)}, fi,) = 

= F(K(r)},{Dfc(r)},{Mfc(r)})- 



ni{r) dr - Ni 



(V) 



E /^[(D.(r)f-A^] dr 

^=1 (V) 



- E / ^ [(M. (r))^ - M^] dr - [ ^r) ( E "^(r) - 1 ) dr 

This functional has the following form: 

/:i({n,(r)}, {M,(r)} , {D,(r)} , {A, (r)}, {u, (r)}, {vl>(r)}, /i,) = 



(39) 



(40) 



where 



i=l 



IJ,N, + ^ / \{r)dr+^ 

(V) iv) 



Ui{r) dr 



+ I \l>(r)(ir 

iv) 



(41) 
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is the part of the functional (139|) . which does not depend in exphcit form on functions Dj(r), 



£1 (K(r)} , {Dfc(r)} , {Mfc(r)}) = F ({n,(r)} , {D,(r)} , {M,(r)}) 



(V) 



E 

i=l 



n,(r) + ^ D,2(r) + ^ M,2(r) + M/(r) n,(r) 



(42) 



is the second part of the functional (139!) . which depends on nj(r), Df (r), Mf'(r) in explicit 
form. The last functional has a form 



(43) 



/:i(K(r)}) = J A(n,(r), Vti,(r)) rfr, 

(V) 

where 'Us(r) denotes all the functions nj(r), -D"(r), M°(r) and 



A(i/,(r), VM,(r)) = --^ 



K^'^ (Vni(r)- Vn,-(r)) + 



+Q!f E [A°(r)n.(r)] ■ V [D;(r)n,(r)]) + 



0=1 
3 



+^!f E [^r(r)n.(r)] ■ V [M;(r)n,(r)] ) 



Q = l 

3 



+•51? E (V [A"(r)^.(r)] ■ V [M;(r)n,(r)]) 



(44) 



4e 



<^ + Ql? (D.(r) • D,(r)) + (M,(r) ■ M,(r)) + 



+4°) (D.(r).M,(r)) 



n,{r)nj{r) + 



1=1 



(E(r) . Di(r)) nj(r) + (H(r) ■ M.(r)) n.(r) + 1 [Aj(r)Dj(r) + v,{r)M^(r)] 



+ E 

2 = 



Tln(^l-/x,-M/(r)a;, 



n r 



The solution of the variational problem for functional fl39|) is equivalent to solution of the 
variation problem for functional fH2l) with conditions ([3]), ([6]), f|T3l) . Thus, the equilibrium 
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distributions of the components and their electric and magnetic moments obey the Lagrange- 
Euler system of equations for functional (P3|) 



dA 



duJr) 



V ■ 



dA 



0, 



(45) 



together with conditions (|3]), ([6]), f[T3l) . 

Using the Lagrange- Euler equation for Us(r) = nj(r) gives the equation: 



E 

i=i 



Kl"^ + (D.(r) ■ D,(r)) + I^^ (M,(r) ■ M,(r)) + Sj"^ (D,(r) ■ M,(r)) 



n,-(r) + 



+ EE (^A"(r) ■ V [Z}J(r)n,.(r)]) + (VMf (r) ■ V [M;(r)n,(r)] ) + 



+^ I^Jf (VDf (r) ■ V [M;(r)n,(r)]) + ^jj^ (VMf (r) • V [DJ(r)n,(r)] ) 



(2) 



+ 



(E(r)-D,(r)) + (H(r)-M,(r))+Tln 



^ rra 



n,- r 



n(r) 



/ii + uji"^ (r) 



ir^) An,(r) + 5^ Qlf (V ■ { A"(r) V {D^{v) n,(r)) }) + 

a 

+1 E '^^^ ( ■ { A" (r) V (M;(r) n,(r)) }) + (V ■ {MfW V {Df{r) n,(r)) }) 



Using the Lagrange- Euler equation for Us(r) = -Df (r) gives the equation: 

J2 Qif ^J(r) + Sif M;(r) n.(r) n,(r) + i?"(r) n.(r) + A.(r)Dr(r)- 



(46) 



^ E 1^;? (^^^(r) • ^ [^"(r)%(r)]) + ^ (Vn,(r) ■ V [M;(r)n,(r)] ) | = (47) 



« V • E ^ [^"(r)^.(r)] + ^ Sj^n^iv) V [M;(r)n,(r 



Similarly, using the Lagrange- Euler equation for Us{r) = Mf (r) gives the equation: 
J2 Rf Mf{v) + - S\f DJir) n,(r) n,(r) + //"(r) n,(r) + z/.(r)Mf (r)- 

-6 E ^ (^^^(^) ■ ^ [^r(r)r^,(r)]) + - (Vn,(r) • V [D^{v)n,ir)]) = (48) 

rn . 

= -g ^ ■ E 1 ^;?^^(^) ^ [^;(r)^.(r)] + - 5(J)n,(r) V [DJ(r)n,(r)] 



The system of equations fH6| l47t |48|) with conditions ([3], O [13]) describes space distri- 
butions of the components and local electric and magnetic moments for the case of full 
thermodynamic equilibrium in the system. 
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C. Partial equilibrium in the system 

It should be noted that in view of colossal times of relaxations for solid state systems, 
the full thermodynamic equilibrium as a rule does not realizable. In the best way, there are 
the particular equilibriums in real solid state structures. 

In particular, for inhomogeneous layered systems with electric and magnetic degrees 
of freedom relaxation times related to electric and magnetic degrees of freedom are short 
whereas times of redistribution of the components have a colossal scales. Thus, the or- 
der parameters related to slow degrees of freedom, are not equilibrium parameters. They 
predetermine by the prehistory of the samples and should be preassigned. 

Suppose that the slow variables in the system are space distributions of the components 
nj(r), whereas distributions of the electric Dj(r) and magnetic Mj(r) moments are the fast 
variables. Then the functions nj(r) should be prescribed, but the functions Dj(r) and Mj(r) 
obey the equations (HH SHD and P El [13]). 



IV. CONCLUSION 

This paper sets as a goal development the unified theoretical approach to both statistical 
thermodynamics and kinetic phenomena in microheterogeneous systems with account their 
most essential real pecuharities: 

1. Presence of some inner degrees of freedom of constituent atoms (such as their atomic 
electric and magnetic moments); 

2. Presence of the competitive short-range and long-range parts of the components in- 
teratomic potentials; 

3. Presence of essential difference of the components atomic sizes. In particular, owing 
to this difference there are no ideal lattices in the real condensed matter. 

It should be noted that such peculiarities as misfits on the phases boundaries, deforma- 
tions fields, and local elastic properties of the matter in microheterogeneous systems are 
caused mainly by the interatomic potentials. Therefore, the correct account of the indicated 
peculiarities of the real systems should be as basis for the adequate model of the real mi- 
croheterogeneous systems, especially for the layered and composite systems with essential 
magnetoelectric interaction. 

There are many works dealing with theoretical research of microheterogeneous systems 
with the magnetoelectric effect. The reviews of some results on the subject my be find in 



papers [15|, Il6|, ll7|, [ISj. Most of the theoretical papers are based on the phenomenological 
Ginzburg-Landau-like free energy functional with additional terms considering lattices dis- 
tortions and misfits on the interphase boundaries [l^, 2^- There is a series of the papers 



by Khachaturyan with co-workers [21|, |22[ based on the phase field approach to composites. 
As a rule the nature of the phenomenological parameters remains hidden. 

It should be noted the very interesting works on the statistical properties and random 
geometry of inhomogeneous systems (23| . 

The present paper contains the other approach to statistical thermodynamics of the mul- 
ticomponent condensed matter. The main peculiarities of this approach consist of following. 
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1. 



Instead of the ideal lattices we introduce the packing condition ([3]). This condition 
allows to take into account essential differences of the components atomic sizes. 



2. We take into account the long-range parts of the interatomic potentials by means of 
the effective fields approximation. Owing to packing condition the GLM as a whole 
falls outside the limits of the effective fields approximation. 

3. We take into account the interatomic potentials depending on inner degrees of freedoms 
such as atomic electric and/or magnetic momentum. The magnetoelectric interactions 
are included in the implicit form. 

4. The Helmholtz free energy functional for the generalized lattice model with account of 
electric and magnetic atomic moments is reduced to GLCH-like form. The connection 
between the interatomic potentials and parameters of the GLCH-like functional is 
established. 

One of the most interesting applications of the GLM is research of layered structures with 
alternation of magnetic and ferroelectric layers. In particular, this approach permits 

• to take into account lattices misfits on the interphases without any auxiliary assump- 
tions; 

• to impart the physical interpretation to the phenomenological models like GLCH; 

• to find the ways for prognosis of the layered structures. 
These results will be presented in the next papers. 
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